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Abstract. We study a simple, solvable model that allows us to investigate effects 
of credit contagion on the default probability of individual firms, in both portfolios of 
firms and on an economy wide scale. While the effect of interactions may be small in 
typical (most probable) scenarios they are magnified, due to feedback, by situations 
of economic stress, which in turn leads to fatter tails in loss distributions of large loan 
portfolios. 
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1. Introduction 

Modelling credit risk in a coherent yet applicable manner is an important yet challenging 
problem. The difficulties arise from the combination of a large, and co-dependent set of 
risk parameters such as default rates, recovery rates, or exposures, which are correlated 
and non-stationary in time. An additional issue is that of credit contagion |T] |2J Hj, 
which examines the role of counter-party risk in credit risk modelling. If a firm is in 
economic distress, or defaults, this will implications for any firm which is economically 
influenced by this given firm, for example, a service provider to it, purchaser of its goods 
or a bank with a credit line to the firm. The direct correlations between firms caused by 
credit contagion lead to further complications in modelling the overall, either portfolio 
or economy wide, level of risk. Jarrow and Yu P introduced a framework of primary 
and secondary firms, the former would default depending on some background stochastic 
process while the latter were affected by a stochastic process and the performance of 
the primary firms. They argued that this was a reasonable level of detail for their 
purposes and it also simplifies matters as there are no feedback loops in the system. 
Secondary firms depend only on primary firms whose performance is independent of 
the secondary firms. Another approach for modelling credit contagion dynamics was 
provided by Giesecke and Weber who used the well known voter process from 
the theory of interacting particle systems, to model interactions between firms. They 
assumed a regular structure for their firms (a regular hyper-cubic lattice) and focussed 
on the equilibrium properties of their model. Davis and Lo jl] considered a model in 
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which defaults occur either directly, or through infection by another defaulted firm, 
with probabilities for direct default or infection taken uniform throughout the system. 
Defaults occurring due to both, endogenous or exogenous causes were not considered in 
their set-up. 

There are a variety of techniques for modelling the correlations between firms' 
default behaviour, which is a major complication in credit risk modelling. The binomial 
expansion technique assumes independence between firms so that the number of defaults 
in a portfolio is described by a binomial distribution. In order to capture the effects 
of correlations a binomial distribution with an "effective" number of firms is assumed 
which is smaller than the actual number in the portfolio, but the weight given to 
each firm scaled so as to keep the mean number of defaults constant, while the 
variance of the overall number of defaults is increased. The relationship between 
the true number of firms and the effective reduced number is a modelling choice that 
depends on the diversity of the firms in terms of sectors, geographic locations or any 
other identifiable trait that would lead to strong correlations in default behaviour. 
JP Morgans' CreditMetrics approach 5J and Credit Suisse First Financial Products 
CreditRisk 4 " |Hj (see [Zj for a detailed comparison between the two) uses the correlations 
in equity values as a surrogate for the correlations in credit quality. The structural 
modelling approach goes back a long way to work by Merton [S] which directly models 
the dynamics of a firm's assets, with default being triggered by the asset value hitting 
some predetermined value (which henceforth we take without loss of generality to be 
zero). Correlations between firms are due to correlations in the dynamics of different 
firms' assets. This approach is very general, as it is relatively transparent to identify 
different driving forces of asset levels and straightforward to include them in the model 
(though the resulting model itself will be non-trivial). However, it suffers from the fact 
that the asset level is not an observable quantity |HJ. On the other hand, the reduced 
form approach gives default rates for a given firm without modelling the underlying 
default process. Correlations are then directly introduced between the default rates. 
There was some discussion in the literature about whether the reduced form model 
could describe the true level of default correlations seen empirically but Yu [TU] seems 
to have answered this question in the affirmative if a suitable structure between the 
default rates is taken into account. 

The approach we take here is a discrete time Markov process (at the microeconomic 
level) where the probability of a default of a given firm in a particular time step depends 
materially on the state of its economic partners at the start of that time step, as 
well as on macro-economic influences. Using techniques developed in the statistical 
mechanics of disordered systems, and recently applied to this specific model in jllj . 
we are able to solve our model exactly, and given our assumptions that we describe 
shortly, this solution takes a particularly simple form despite the fact that in principle we 
have feedback correlations, non-equilibrium dynamics and in principle non-Markovian 
behaviour at the macroscopic (economy/portfolio wide) level. We note that it is possible 
to frame our model in either the structural approach or the reduced form approach, 
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depending on requirements and taste, although the interpretation of the variables in 
the two approaches will of course be different. We find that the correlations introduced 
through credit contagion lead to large increases in default rates in times of economic 
stress, above and beyond those introduced by simple macro-economic dependencies. 
This has strong implications for portfolio risk management. 

2. The microeconomic framework 

We will analyse an economy of iV firms in the large N limit. Generally, we focus on the 
characteristic changes in the economy due to interactions between firms, which will be 
described in a probabilistic manner. 

As mentioned in the introduction we take a discrete time approach. For clarity 
we restrict our discussion to a one year time frame split into twelve steps; this is not 
essential, but parameters may need rescaling depending on the set-up. We use a binary 
indicator variable n iyt to denote whether firm % is solvent at time t {n^t = 0) or has 
defaulted (n^t = 1). The default process is a function of an underlying stochastic 
process for each firm in terms of a "wealth" variable W^t, where we assume default if 
the wealth drops below zero. We shall assume that recovery from default over the time 
horizon of a year is not possible, so that the defaulted state is absorbing. As a function 
of the wealth, therefore, the indicator variables evolve according to 



where 0(. . .) is the Heavyside function. 

A dynamic model for the indicator variables is obtained from by specifying the 
underlying stochastic process for the wealth variables W^t- We shall take it to be of the 
form 



Here $j denotes an initial wealth of firm i at the beginning of the risk horizon, and 
quantifies the material impact on the wealth of firm i that would be caused by a default 
of firm j. This may or may not be a reduction in wealth, depending on whether j has 
a cooperative (Jjj > 0) or a competitive (Jy < 0) economic relation with i. 

We shall assume that the fluctuating contributions 77^ to (j2J) are zero-mean 
Gaussians. There is still some degree of flexibility concerning the decomposition of 
the rjij into contributions that are intrinsic to the firm and extrinsic contributions. The 
latter describe the influence of economy-wide fluctuations or fluctuations pertaining 
to different economic sectors, depending on the level of detail required. We restrict 
ourselves to a minimal model containing a single macro-economic factor (assumed to be 
constant over a risk horizon of a year), and individual fluctuations for each firm, 




(1) 



N 
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where a, sets the scale of the individual fluctuations, and the {Ci,t} are taken to be 
independent jV(0, 1) Gaussians; finally, the parameters pi quantify the correlations of the 
rji^t created via the coupling to economy-wide fluctuations i] , also taken to be Af(0, 1). 

Up to this point the wealth dynamics does not contain an endogenous drift. If 
predictions are required over longer time periods then it may also be pertinent to 
introduce such a drift, e.g. by using a time- dependent $j for example, # i)t = $j(0)e Zi *, 
where z± denotes an intrinsic growth rate of the average wealth of firm i (with z; t > 
for a firm making profits and Z{ < for a firm making losses). However, for the current 
purposes of examining default rates over the medium term and especially focussing on 
the behaviour on the tails, this adjustment does not lead to significant changes in our 
overall conclusion. 

The model, as formulated above, clearly takes a structural point of view on 
the problem of credit contagion. However, we note that the dynamics (jlj of the 
indicator variables is clearly independent of the scale of the wealth variables W^- By 
appropriately rescaling the initial wealths di and the impact parameters Jy we can 
thus assume a unit-scale crj = 1 for the noise variables Q. Interestingly, this simple 
rescaling, which leaves the dynamics of the system unaffected, amounts to changing to 
a reduced-form interpretation of the dynamics. 

To see this, note from (J2J) that the event W^ t < is equivalent to f] itt > 
$i~Y,jLi Jij n j,t- With o"j = 1, we see that this occurs with probability &{J2j Jij n j,t — 'di) 
where $(■) is the cumulative normal distribution. From a reduced form point of view this 
is just the intensity of default of firm i at time step t (in a given economic environment 
specified by the set of firms defaulted at time t). This allows us to re- interpret the 
(rescaled) initial wealth and impact variables di and in terms of the bare default 
probabilities [TT1 [TJJ [T3] . I.e., if company i has an expected default probability of Pi in 
a given time unit (e.g. one month in the present set-up) as predicted from tables from 
ratings agencies, then di = — Similarly, the expected default probability Pi\j of 
firm i, given that only firm j has defaulted leads to the value = $ _1 (pj|j) — $~ 1 (]?;). 

In determining the model parameters by the method suggested above we are 
splitting our default probability into terms that come from credit contagion and other 
terms such as the bare default probability that come from historical data. It could 
fairly be argued that the historical data already incorporate the credit contagion terms 
and thus we are double counting. As we will see later in numerical simulations, the 
credit contagion terms make very little difference to average behaviour and thus making 
estimates based on average historical data is still a reasonable approach. 

In choosing the variable pi we follow the prescription given by BASEL II ^1] which 

sets 

i _ p -50PDi / i _ -50PD* \ 

a = o- 12 -T^r + °' 24 \ l - T^r) K °' 12 i 1 + e "° PD ') (4) 

where PDj gives the probability of default of firm i over one year, ignoring credit 
contagion effects. With p^ = $(— as the monthly default probability, we have 
PD,- « 12$(-#A 
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We still have to specify the form for the economic interactions. We adopt here a 
probabilistic approach, so take them to be random quantities of the form 



C + J?* 



(5) 



Here, the G {0, 1} detail the network (presence or absence) of interactions between 
different firms and we choose these to be randomly fixed according to 

p ( c ij) = ^f S ai,i + i 1 ~ ^) s cafi > * < 3 , Cji = . (6) 

We assume that the average connectivity c of each firm is large in the limit of a large 
economy; this will allow the influence of partner firms to be described by the central 
limit theorem and the law of large numbers. Concerning the values of the (non-zero) 
impact parameters, we parametrise them as shown, with x^ assumed to be zero-mean, 
unit-variance random variables, with finite moments, and independent in pairs, 

Xi] = , x\a — 1 , XijXji = a , XijXki = otherwise . (7) 

The parameters Jq and J determine mean and variance of the interaction strengths; the 
scaling of mean and variance with c and sfc respectively in (jSJ) is necessary to allow 
a meaningful large c limit to be taken. Taking Jo > would encode the fact that on 
average firms have a synergy with their economic partners. 

At first sight, specifying the J^- appears to introduce a vast number of parameters 
into our model, but in fact only the first two moments of the distribution of interaction 
strengths are sufficient to determine the macroscopic behaviour of the system, and so 
the model space is not too large. 

Let us now turn to the capital required to be held against credit risk. In the BASEL 
II document ^1] the capital requirement for a unit-size loan given to firm i is 



Kt = LGD, 



$ / yp-$- 1 (0.999) + $- 1 (PD t ) \ _ pD? 



Mi . (8) 



The first factor, the loss given default LGD, of firm i, is related to the average fraction 
of a loan that can be recovered despite default. The last factor, Mj, is related to the 
maturity (long dated loans are inherently riskier). Adjustments related to liquidity (low 
liquidity loans are risker) and concentration (fewer, larger loans give a greater variance 
in returns for given expected return) are occasionally also included in this factor - 
concentration-adjustments, in fact, are a means to account for reduced granularity in a 
credit portfolio resulting from the possibility of credit contagion. 

The factor inside square brackets in (JBJ) is entirely related to the loss-frequency 
distribution. The first term is the value of the loss frequency not exceeded with 
probability q = 0.999 under fluctuating macro-economic conditions, with pi describing 
the dependence of the firm's loss-frequency on the macro-economic factor. The second 
term is the average loss frequency. The value of the confidence level q is in principle 
arbitrary, but is related to the target rating of the bank. The risk weighted asset is 
then found by further multiplying by terms such as the exposure at default (i.e. size of 
the loan). Thus the capital required for firm i can be viewed as the loss at the 99.9th 
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percentile level of stress, in excess of the expected loss, multiplied by a conversion 
factor. From this structure it is clear that a key ingredient for the capital adequacy 
requirements is a good model of credit risk that works well into the tail of the loss 
frequency distribution. 

Returning to our description of default dynamics, let us first focus on the case of 
independent firms, with = Vi,j, and consider a single epoch for our model with 
fluctuating forces given by Q at given macro-economic condition r/ . The probability of 
a default of firm i with average unconditional monthly default probability pi occurring 
during the epoch t — > t + 1 in our model is given by 

(^,mKt = 0) = $ \ - — _ I (9) 

Since the probability of default is increasing with 770, we can find the probability of 
default not exceeded at e.g. the 99.9 percent confidence level; it is given by setting 
7/0 = (0.999) in the above equation (recall r/ is distributed as a zero-mean, unit- 
variance Gaussian random number). As above, the excess capital required is the loss at 
the 99.9th percent level minus the expected loss (multiplied by a risk factor). However, 
when we consider the case of an interacting economy with non-zero J^, we find that in 
fact 



V vi - Pi + J m ) 

m t = jjH n M (11) 

3 

is the fraction of firms within the economy that have defaulted up to time t; we also 
expressed the expected monthly default rate pi in terms of a 'rescaled initial wealth' $j, 
$-i(p.) = 

Thus we find that our formulation is very similar to that used in BASEL II. However, 
we directly take account of the correlations in defaults caused by credit contagion. 
This introduces two extra parameters into the model but it does markedly change the 
behaviour in the tails of the loss frequency distribution, and thereby in the tails of the 
loss distribution itself. Correlation between firms is essentially a dynamic phenomenon 
- if there is no dynamics, there is no way for one's firms' performance to influence 
the viability of any other firm. Thus rather than considering firms to be independent 
over a single epoch which lasts the entire period of any loan, we split the overall time 
(e.g. one year) into smaller units (e.g. one month) and let the firms evolve over these 
smaller time units with the default probability adjusted (since the default event in 12 
monthly epochs is compounded 12 times as opposed to a single epoch). A firm may 
default at any point, but will then influence its partner firms for the remainder of the 
time horizon. The complexity of the theory is merely linear in time, thus it is not a 
great computational burden to choose this approach. 

Following the approach described in [TJ] it is possible to solve the model in a 
stochastic manner. Credit contagion within this model is encoded at each time by a 



where 
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single number, the fraction of firms that have defaulted thus far, which evolves according 
to 

J m t + Jp~(¥) 7 l0 - 

>/i - + ^ )K 

where (n t {"d)) denotes the time-dependent monthly default rate of firms with ^ pa as 
influenced by interactions with the economy (see Eq (16) of and the larger angled 
brackets with subscript d denote an average over the bare monthly probabilities of 
default for the ensemble of firms, or equivalently over the distribution p{"&) of their 
rescaled initial wealth parameters In (|12|) the Basel II recommendation which 
links correlations to macro-economic factors with (unconditional) default probabilities, 
Pi = p(pi) — > p(i?i), via (JH) is formally taken into account. 

Note that credit contagion affects the dynamics of defaults only via two parameters, 
Jo and J, which characterise the mean and variance of the impact parameter 
distribution. Note also that the parameter a which quantifies forward-backward 
correlation of mutual impacts according to does not appear in the final formulation, 
nor are there any memory-effects in the dynamics, as would normally be expected for 
systems of this type. The reason for this simplifying feature is in the fact that the 
defaulted state is taken to be absorbing over the risk horizon of one year. 

3. Results 

We now turn to presenting a few key results of our analysis. Our results concerning 
default dynamics and loss distributions are obtained for an economy in which the 
parameters #j determining unconditional monthly default probabilities Pi according to 
i?j = — $ _1 (pj), are normally distributed with mean i9 — 2.75, and variance a\ = 0.1 
so that typical monthly bare default probabilities are in the 10~ 5 — 10~ 3 range. The 
couplings to the macro-economic factor are chosen to depend on the expected default 
probabilities according to the Basel II prescription (jlj). 

In Fig. 1 we we show that renormalisation (with respect to credit contagion) 
makes little difference to the typical default dynamics. The evolution of the fraction of 
defaulted firms in interacting economies differs hardly from that of the non-interacting 
economy with J^- = <=>■ ( J , J) = (0, 0). 

In marked contrast to this, the tails of loss-frequency and loss distributions are 
strongly affected by the presence of interactions in the system, as shown in Figs. 2 and 
3. We note that the tails of the loss-frequency distribution and the loss distribution 
are more pronounced than in our previous study This solely due to the fact that 
in the present paper we followed the Basel II suggestion that relates the coupling of a 
company to macro-economic factors with its default probability via (jlj). 

In computing the economy-wide losses per node for given macro-economic condition 

Vo, 

L{vo) = -rz XXi2^ (13) 



m t+ i = m t + 



<n t (tf)))$ 
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we assume that the ii are randomly sampled from the loss distribution for node i, taken 
to be independent of the stochastic evolution, but correlated with the bare monthly 
default probability. In the large N limit this gives 

L( Vo )= lim i£n 12 (^)4= f dtipWinuWW) (14) 

N— >oo iv J 
i 

by the law of large numbers, where J = is the mean of the loss distribution 

for a node with default probability Pd{$)- As an example we consider an economy 
where average losses are inversely proportional to the unconditional default probabilities 

Pi=Pd{$i) = 

£(■&) = £ ±— (15) 

with a parameter e preventing divergence as Pi —>■ 0. In this way, the contribution 
to the total losses will be approximately uniform over the bands of different default 
probabilities. The distribution of the economy-wide losses per node is driven by the 
distribution of the macro-economic factor, and is computed analytically as shown in 
|llj . A typical result is shown in Fig. 3, for which we chose the scale l§ — 1 and the 
regularizer e = 0.005. Once more economic interactions are seen to strongly affect the 
tail of the loss distribution at large losses, which is due to the possibility of avalanches 
of loss events in times of extreme economic stress. 

Note that we have been dealing here with "synthetic" parameter distributions for 
averages of loss distributions, as well as bare monthly failure probabilities. These could 
be replaced by realistic ones without affecting the general set-up. We have not looked 
specifically at finite size effects here. In it was shown that they are fairly small. 
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Figure 2. Distribution of the fraction of defaulted companies for (Jo, J) = (0,0) 
(bottom) and (J , J) — (1, 1) (top). 
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Figure 3. Distribution of losses for the non-interacting system with (Jo, J) = (0,0) 
(bottom) and the interacting economy with (Jo, J) = (1, 1) (top). 

4. Conclusion 

In this paper we have looked to incorporate the risk due to credit contagion into 
the internal ratings based approach discussed in BASEL II. While the mathematical 
subtleties are discussed in full detail elsewhere fTTj . essentially the large number of 
neighbours assumed for firms means that the law of large numbers and central limit 
theorems apply to the interactions, meaning that our theory requires only two more 
parameters than the BASEL II approach. In terms of risk, one of the striking results 
is that while the effect of interactions is relatively weak in typical economic scenarios, 
it is pronounced in times of large economic stress, which leads to a significant fattening 
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of the tails of the portfolio loss distribution. This has implications on the fitting of loss 
distributions to historical data, where care must be taken not only to fit the average 
behaviour but also to take care with the more extreme events. 

This touches the issue of model calibration discussed in greater detail in [T3] . 
We note that our model requires bare default probabilities and conditional default 
probabilities as inputs. Historical data, however only contain interaction-renormalised 
default probabilities, and thus the problem arises of how to disentangle the two effects. 
Concerning typical behaviour, Fig. 1 shows that the effect of interactions is fairly small, 
and interaction-renormalised default probabilities can, to a first approximation within 
this model, be taken as substitutes for the bare ones. Concerning conditional default 
probabilities, these would have to be obtained from refined rating procedures; see |13j . 
Interestingly, however, only the low order statistics of these are needed to describe the 
collective dynamics of the system. Their effect manifests itself only in situations of 
economic stress, generating fat tails in portfolio loss distributions. 

The model we have proposed is relatively simple in two important respects. Firstly, 
we do not take into account credit quality migration but have just two states for our 
firms, solvent or defaulted. The model could be extended to allow for more states for 
each firm, although the full complexity of non-Markovian dynamics would resurface in 
an attempt to take credit quality migration along these lines into account. Secondly, 
the firms and their environment are rather homogeneous, which in practical situations is 
of course an approximation. This approximation has been made for convenience rather 
than out of necessity; the techniques described in [TT] can be adapted so as to treat 
situations with more heterogeneity in local environments. We intend to work on some 
of these possible model generalisations in the future. 

One advantage of our simple model is that it is exactly solvable and the solution 
itself is not overly involved theoretically or computationally, and we only need to 
introduce two extra parameters to quantify the effect of economic interactions - 
compared to the BASEL II approach, which ignores credit contagion altogether. 
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